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Abstract 

The elementary local and global influence of geodesic field line curvature on radial dispersion 
of zonal modes in magnetised plasmas is analysed with a primitive drift wave turbulence model. 
A net radial geodesic forcing of zonal flows and geodesic acoustic modes can not be expected in 
any closed toroidal magnetic confinement configuration, since the flux surface average of geodesic 
curvature identically vanishes. Radial motion of poloidally elongated zonal jets may occur in the 
presence of geodesic acoustic mode activity. Phenomenologically a radial propagation of zonal 
modes shows some characteristics of a classical analogon to second sound in quantum condensates. 
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I. INTRODUCTION 



Zonal modes play an important role in turbulent self organisation of quasi two- 
dimensional fluids, such as magnetised plasmas or in geophysical, planetary and stellar fluid 
dynamics [1-3]. Zonal flows can be regarded as a low- frequency spectral condensate phase 
of the turbulence [4, 5] featuring the highest flow symmetry allowed by the (generically in- 
homogcncous) system: atmospherical and oceanic zonal jets are latitudinally extended flow 
structures, and zonal modes in toroidal fusion plasmas appear as average flows on magnetic 
flux surfaces. 

Geophysical zonal flows can experience temporal variations and oscillations caused by 
seasonal or topographical influences, and may migrate meridionally. Similarly, plasma 
zonal flows in the inhomogeneous toroidal magnetic fleld are modulated by flnite frequency 
geodesic acoustic mode (GAM) oscillations [6]. The possibility of radial propagation of zonal 
flows and GAMs in fusion plasmas has recently attracted interest, motivated by indicating 
experimental observations [7-9]. 

Previous studies on radial propagation included analytical theory of GAMs and zonal 
flows [10-12] and linear simulations of GAMs [13-15]. Finite Larmor radius (FLR) effects 
in the presence of a background temperature gradient have been identifled as one potential 
mechanism for radial propagation of zonal modes [13-17]. Further, it had been postulated 
that details of the magnetic field geometry, in particular an up-down asymmetry, could cause 
radial GAM migration [18, 19]. 

In the following the particular influence of geodesic curvature and up-down asymmetry 
on the radial dispersion of zonal modes is analysed. A basic geometric property has to be 
recollected in this context: the flux surface (zonal) average of geodesic curvature is zero in 
any toroidal configuration. Real and artificial geodesic forcing effects on zonal modes are 
discussed, including the perils of using inadequately simple geometry models that violate 
this property. 

The computational analysis shall here be based on a self-consistent nonlinear model in- 
cluding both the drift wave turbulence ffuctuation component and the zonal mode conden- 
sate. Radial propagation can best be studied when radial boundary conditions are eliminated 
by using a periodic domain. This however excludes the utilization of 3-D ffux-tube models, 
which rely on finite magnetic shear to avoid unphysical flute modes but thereby introduce 
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radial secularity. 

The most primitive model for self-sustained drift wave turbulence in magnetised plasmas, 
that also allows radially periodic boundary conditions, is the two-dimensional Hasegawa- 
Wakatani (HW) set of equations [20]. This standard model is here extended to specifically 
study geodesic curvature effects. 

The paper is organised as follows: in section II the curvature modifications to the HW 
model are introduced. The linear local radial dispersion relation of geodesically forced zonal 
HW modes is discussed in section HI. Numerical details regarding the nonhnear simulations 
are presented in Chapter IV. Computational results on local geodesic forcing of turbulence 
generated zonal flows are reported in section V. The essential flux-surface property of zero 
average geodesic curvature and pitfalls concerning simplified geometry models arc reviewed 
in section VI: violation of this constraint results in unphysical radial propagation of zonal 
flows and GAMs. Poloidal variation of curvature is introduced in section VII, and effects on 
apparent radial propagation of zonal modes is discussed in section VIII. Finally, in section 
IX, a striking analogy to the phenomenon of "second sound" in quantum condensates is 
presented. 

II. TURBULENCE-FLOW MODEL 

First, consider the HW model for drift wave turbulence [20]. It accounts for nonlinear 
instability driven by a gradient Vno(x) in plasma density and resistive parallel coupling 
between fluctuations of density n and electrostatic potential (j). The resulting turbulent 
state oi E X B vortices in the (x, y) drift plane perpendicular to the magnetic field can form 
low- frequency ky = zonal flow structures with a finite wave number k^. The HW model 
is in the following employed with modifications for zonal flow corrections on the dissipative 
coupling [22] and including magnetic field line curvature terms [23]: 



Standard drift normalisation is applied [23]. The electrostatic potential is obtained from 
the vorticity Q by 



dtn + [(p,a\ =d{(p 

dtn + [0, n] = d{(f) 



h) — K{n), 

- gndy(f) + K{(l)-n). 



(1) 
(2) 



(3) 
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In the dissipative coupling term the zonal (y) average is subtracted from density and 
potential fluctuations [22]: 

4) = (j)—{4>), h = n — {n). (4) 

In toroidal plasma geometry the x coordinate locally represents the minor radial, and y the 
poloidal direction. 

The original HW model (with k = 0) is determined by the dissipative coupling param- 
eter d (proportional to /c^) and the density gradient scale length Qn- The hydrodynamic 
limit of the Euler equations is recovered for c? = and Qn — 0, while the adiabatic limit 
1 asymptotically corresponds to the Hasegawa-Mima-Charney-Obukhov equation. The 
general properties of the HW model have been extensively discussed elsewhere (e.g. in 
Refs. [24-33]). 

The extended model here further includes held line curvature effects. Normal and geodesic 
components of the magnetic curvature /C = Vlni? enter the compressional effect on vortices 
due to magnetic field inhomogeneity by 

K — Kiydy + Kxdx (5) 

where the curvature components in toroidal geometry are a function of the poloidal angle 
0. For a circular torus Ky = cbCOs(^) and = CBsin(^) when ^ = is defined at the 
outboard midplane. Pure interchange drive is acting for ^ = and pure geodesic effects 
for 9 — ±7r/2, and mixed effects are achieved for intermediate angles. Slab turbulence is 
recovered for cb = 0. 

In the following, at first the local influence of geodesic curvature on zonal modes in drift 
wave turbulence is discussed by choosing a constant angle 9 for the whole y-domain. Later, 
a periodic dependence of the form 9{ii) — 27: {y — Ly/2)/Ly will be introduced to globally 
map the magnetic curvature inhomogeneities on a flux surface into the simulation domain. 



III. LINEAR ANALYSIS: LOCAL GEODESIC FORCING OF ZONAL MODES 

Fourier analysis with an ansatz n, ^ ~ exp(ik • x + iuit) and linearisation of eqs. (1) and 

(2) for ky <^ kx and o? = (zonal modes) in the local approximation provides 

ukl(f)k = i^xKrik (6) 
ujUk = KxKih - rik) (7) 
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For global simulations covering the whole flux surface (in an annulus or flux-tube repre- 
sentation) the averages {Kxnk)y and {Kx(f)k)y oi the geodesic curvature operator acting on the 
fluctuations would have to be taken into account in the linear analysis (see section VII): side- 
band modes of the density then will result in the well-known geodesic acoustic oscillations 
at the GAM frequency. 

In the present local discussion in this section, {Kx)y in itself is assumed to be nonzero 
and constant, and global sideband modes are neglected. This artiflcial scenario however can 
apply to zonal jet modes with nonzero but small ky <^ kx (rather than zonal flows with 
exactly ky = 0), and it corresponds to artiflcial situations where a flnite flux-surface average 
geodesic curvature results from an improper geometry model. 

The resulting local zonal wave equation has a frequency 

uj ^ Cxkx{n ± \) (8) 

which is nonzero for flnite average geodesic curvature, where = —{kx/^tt) and A = 
{tt / kx) ^/^+~kl. The flnite kx zonal modes are actually nonlinearly driven by turbulent 
Reynolds stress, which is eliminated in the zonal dispersion relation eqs. (6) and (7) by 
linearisation. 

The dispersion relation results in a phase velocity 

Vph^ ^ ^ Cx{7r ± X). (9) 

For kx <^ 1 follows \^2T^/kx so that u ^ 27rCx is asymptotically independent of the wave 
number. Radial modulation of zonal modes is in this long wavelength limit obtained by a 
local phase velocity 

Vph ~ CxX. (10) 

The group velocity Vgr — dui/dkx ~ ttCx is for kx <^ 1 much smaller than the phase 
velocity. The following numerical simulations will show that zonal flows and GAMs mostly 
appear with one dominant radial mode number propagating with Vph. For single modes the 
concept of a wave packet travelling with a group velocity is not relevant. Thus a Poynting 
approach on GAM propagation [18] is pointless. 

In the following the simple analytical relation (10) is numerically tested locally, and the 
effect of an artiflcial non-zero average geodesic curvature Cx on radial propagation of zonal 
modes is demonstrated and visualized. 
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IV. NONLINEAR NUMERICAL SOLUTION 



Equations (1) and (2) are numerically solved in time with an explicit 3rd order Karni- 
adakis scheme [39], and the Poisson brackets [a, 6] = {dxa){dyb) — {dya){dj)) are evaluated 
with the energy and enstrophy conserving Arakawa method [40]. The numerical method is 
equivalent to the one in Refs. [23, 32, 41]. Hyperviscuous operators z/^V^, with = — 2-10~^, 
are added for numerical stability to the right hand side of both equations (1) and (2), acting 
on Q and n, respectively. 

Eq. (3) is for a double periodic domain efficiently solved in k space by evaluation of 
0fe = —Qk/k±. Threaded FFTW3 hbraries are employed for the 2-D forward and backward 
Fourier transforms [42]. The use of openMP multi-threaded FFTW3 routines on an 8-core 
workstation turned out to be parallel efficient up to 4 threads only for the specified grid size 
due to overhead costs. 

The equations are here discretised on a doubly periodic grid with various (in general not 
quadratic) box dimensions. Spatial scales are in units of the drift scale Ps — ^/T^Mi/ i^B) 
where Te is the electron temperature, Mj is the ion mass, and B is the magnetic field strength. 
The gradient length scale Qn — L±/Ln = 1 is fixed by specifying L± — — iVlnriol"^. 
This normalises the time scale to t ^ t Cg/ Ln- A stable time step for nominal parameters 
on a 512 x 512 grid with = Ly = 128ps scale was found to be At = 0.0025. The 
computations are initialized with a random bath of quasi-turbulent density fluctuations, 
run into a saturated turbulent state (typically obtained for t > 500) and are continued for 
statistical analysis. 

V. LOCAL RESULTS: RADIAL PROPAGATION OF ZONAL MODES 

The radial zonal flow structure with respect to local normal and geodesic curvature effects 
is in the following numerically analysed by varying the poloidal angle 6 through different 
simulations, while 9 takes a constant value on the respective {x, y) domains. 

Nominal values for the dissipative coupling and curvature parameters used here are d — 
0.5 and cb — 0.05. The flows are dominantly zonal for these parameters, and are blocked for 
lower forcing [d ^ 0.1 and without interchange drive). The radial mode number I of zonal 
flows increases with stronger forcing. 
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FIG. 1: Radial wave number kx (in units of the drift scale ps) spectra of the potential Bold 
line: mixed curvature effect (cb = 0.05, 9 = tt/8); thin dashed line: slab case (cb = Oj. 

Radial wave number spectra I^KA;^;) for the d = 0.5 flow dominated scenario are shown 
in Fig. 1: For the standard HW slab case {cb = 0, thin dashed line) the radial wave number 
features a distinct peak at ~ 0.39 in units of Ps- This corresponds to an / = 8 radial 
mode of the {m — 0, n — 0) zonal flow with radial wavelength A = li^/kx ~ 16ps. Further 
contributions from I — 6 and I — 10 modes are an order of magnitude smaller. The mixed 
curvature case {6 — tt/S, bold line) shows a strong maximum for I — 10 and a broader 
spectrum with similar amplitudes for I = 4, 5, 7, and 8 modes. 

This radial structure of the zonal flow modes is also clearly visible in 2-D {t,x) plots for 
cases with strong forcing. Fig. 2 shows the time evolution of the zonally averaged potential 
for various drive parameters and poloidal angles. Here the spatial resolution is reduced to 
256 X 256 grid points for the same physical domain size of Ly — 128ps as before. Runs are 
now taken to t — 10^. 

As predicted above, the zonal modes propagate radially with a velocity proportional to 
the geodesic curvature parameter and to the wavelength A = L^/l. The phase velocity 
measured in the nonlinear simulation is in very good agreement with the estimate obtained 
from linear theory. For case (e) with Cx — —CB/i^n) ^ —0.008 and X/ps — 128/8 — 16 we 
expect Vph ~ CxX ~ 0.064 from equation (10), while the simulation gives Vsim — Ax/ At = 
128/1860 fa 0.069. 
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FIG. 2: Zonal potential {(f)){t,x) for various cases of forcing and curvature parameters: 
(a) blocked Bow for low forcing (d = 0.1, cb = 0); (b) zonal flow for high forcing (d = 0.5, 
cb = 0.05, 6 = 0); (c) mixed curvature case (d = 0.5, cb = 0.05, 9 = ir/8); (d) same but reversed 
field (d = 0.5, cb = 0.05, 6 = -n/8); (e) pure geodesic drive (d = 0.5, cb = 0.05, 9 = 7r/2). 



At this point it should again be stressed that this local result as a test for eq. (10) does 
not apply to actual toroidal ky = zonal flows. It will however be relevant for the following 
discussion on geometrical artefacts and on the interpretation of spatio-temporal features 
identified as radially propagating zonal jet structures. 
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VI. FLUX-SURFACE PROPERTY OF GEODESIC CURVATURE 



In a torus the "upper" and "lower" regions locally correspond to positive and negative c^, 
respectively, if the poloidal coordinate 9 is chosen as here to be aligned "upwards" . Reversal 
of the magnetic field direction will locally also reverse the sign of the geodesic term through 
Cx{—d) = —Cx{0). For a circular torus Hx = CBsm{9). 

The poloidally local effects of radial geodesic forcing globally balance, if positive and 
negative contributions from geodesic curvature cancel by integration over a flux surface. 
This is indeed the case for any low beta toroidal MHD equilibrium, including up-down 
asymmetric tokamaks or 3-D stellar ator configurations. The flux surface average of the 
geodesic curvature always identically vanishes [34]. 

This is illustrated in Fig. 3 (a) where the geodesic curvature Kx{d) as a function of the 
poloidal angle is shown for the example of a real up-down asymmetric ASDEX Upgrade 
equilibrium in a lower single-null configuration, calculated with the VMEC equilibrium code 
[35] at a radial position p — V/Vg — 0.95 in volume coordinates. The maximum of geodesic 
curvature is reduced in the vicinity of the X-point location 9 — +7r/2, but the average over 9 
remains zero as expected. This peak reduction affects the GAM and zonal flow amplitudes 
(and consequently turbulent transport) through the same mechanisms as by the up-down 
symmetric effect of ffux-surface elongation [36-38]. 

In no toroidal configuration a net radial propagation of (m = 0, n = 0) zonal modes by 
geodesic forcing alone is therefore possible. 

Approximate analytical toroidal geometry models (e.g. when introducing up-down asym- 
metry effects) have to account for this cancellation, otherwise unphysical radial propagation 
may appear in simulations. This is illustrated in Fig. 3 (b): here k,x{9) (bold line) has been 
calculated from the simple up-down asymmetric model geometry from Ref. [19], where flux 
surfaces are vertically shifted. It can be clearly realised that this corresponds to a shift 
Kx{9) ^ Co + cb sin 9 with the result of a nonzero average geodesic curvature. This unphys- 
ical finite average acts on GAMs and zonal fiows exactly like the local model discussed in 
Sections III and IV and results in an artifical radial propagation. Other simplified analyt- 
ical geometry models in general also do not conserve the zero average geodesic curvature 
property. 

Most results concerning the overall turbulence and fiow levels and other characteristics 
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FIG. 3: (a) Geodesic curvature Kx{0) in a circular torus (sinusoidal thin curve) and in a realistic 
ASDEX Upgrade equilibrium (bold curve). The average always is identically zero, 
(b) Geodesic curvature Kx{0) resulting from an inappropriately simple X-point geometry model 
(bold curve), which leads to an unphysical hnite average. 

from nonlinear simulations arc usually not or only to a very small degree affected by this 
particular inconsistency. Simple X-point models for example well allow to study an influ- 
ence of local magnetic shearing on turbulence and transport. For any discussion of radial 
propagation of GAMs and zonal flows the use of analytical model geometries has to either 
explicitly account for a zero average or is otherwise inappropriate. Then numerical solutions 
from codes that solve actual MHD equilibria are required. 
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VII. 2-D GEODESIC ACOUSTIC MODES 



Next, a poloidal variation of the geodesic curvature within the 2-D HW computational 
domain is taken into account. The extension Ly is elongated up to Ly/L^ — 16 in order 
to cover a complete poloidal circumference. The curvature terms Ky = cb cos{6) and = 
CBsin(^) are now periodically varying along the y coordinate by 9{y) = 2n{y — Ly/2)/Ly. 
The average along y can so be interpreted as the flux surface average, while still any parallel 
(3-D) dynamics is neglected. 

The hnear spectral zonal HW equations now become 

^Ki'Pk) = (f^xKrik) (11) 
u{nk) = {nxkx{.(t>k - rik)) (12) 

where the brackets denote the average along y. Using the identity sin^ — \ — cos^ 9 — 
|[1 — cos(2^)] and defining the zonal flow velocity by = ikx{(l>k), the dispersion relation 
gives cu'^Uk = \c\uk. The result is a periodic oscillation of the zonal flow with the GAM 
frequency uj — CBj\f2 = uqam- 

The present 2-D HW model lacks any parallel effects of field line connection and wave 
dynamics. In particular, the absence of coupling from the density sideband to parallel 
dissipative sound and Alfven dynamics overestimates the GAM amplitude, compared to 
3-D fiux-tube simulations [43]. On the other hand, the density sideband is also in 2-D 
nonlinearly coupled to the density cascade and its respective small-scale dissipation. The 
geodesic transfer thus constitutes a sink mechanism for zonal fiow energy, and the 2-D 
toroidal zonal fiow amplitude is less than in the slab or local case (and also less than in 3-D 
toroidal simulations). 

The results are clearly quantitatively different from more complete 3-D (drift-Alfven or 
gyrofiuid) models, but the emphasis here is on principle understanding of basic curvature 
coupling mechanisms which do not specifically rely on parallel dynamics. 

VIII. SPATIOTEMPORAL ZONAL FLOW PATTERN: APPARENT RADIAL 
MOTION 

The poloidally extended model for k — K{y) is now solved numerically with simulation 
parameters d — 0.5, — 1, Ux — 128, Uy — 2048, — 32ps, Ly = ^Vlps, — —b ■ 10~^, 
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FIG. 4: Pattern of zonal Bow velocity {vy){x,t) for cb = 0.5 (top) and cb = 0.25 (bottom). 
Intermittent diagonal streaks bear semblance of radial "propagation" for high-frequency GAMs at 
Cb = 0.5 but are absent for lower frequencies. 

At = 10~^, and with a curvature parameter cb = 0.5. The zonal flow velocity {vy){x,t) is 
shown during the initial saturated phase of the computation (100 < t < 900) in Fig. 4. 

The plot for cb = 0.5 features an irregular checkerboard like pattern: quasi periodic 
oscillations appear both in space and time as the GAM modulation reverses the zonal flow. 
The radial wave number of the pattern is in the order ^ 1 and the frequency in the 
order of w ~ coqam- When the curvature parameter is reduced to cb = 0.25 (with otherwise 
identical simulation parameters) the frequency scales proportionally but the modulation 
appears weaker and the zonal flow direction at a specific radial location x rarely reverses 
with time. 

Another feature of the pattern for cb = 0.5 immediately strikes the eye: some zonal flow 
structures appear connected diagonally and give the semblance of moving radially in time 
with a velocity f"^^ = u/k^. The connections first appear randomly with no prefered radial 
direction. This in-out symmetry can be spontaneously broken and an effective diagonal 
pattern may emerge transiently (as can be seen for around t ~ 100 and t ^ 500). Occur- 
rence and directionality of such symmetry breaking appear to be random and to depend on 
simulation parameters and resolution. 

The resulting apparent radial motion however is illusionary as the actual spatio-temporal 
zonal checkerboard-like structures do not move. This appearence may mislead to the con- 
clusion of GAMs propagating radially with v^^^ ~ ^gam/^x- But as it turns out it is not the 
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FIG. 5: "Zonal jets": ky <^ 1 vortices strained out by zonal Bow shearing (in particular around 
X ^ 50). Vorticity i^{x,y) is shown (negative values: blue; positive values: red; in arbitrary units) 
on a part of the computational domain. 

GAM modulated zonal flows (or the GAMs) which are moving, but rather ky 1 vortices. 

Vortex structures can be detected in a wide spatial range, from circular eddies in the size 
of ky ^ kx ^ 1 to ky = zonal flows. In an intermediate range poloidally elongated vortices 
with ky <^ 1 are visible. These "zonal jet" structures appear strained out by the zonal flow 
shear, like shown in Fig. 5. 

Zonal jets are, in contrast to ky = zonal flows, susceptible to geodesic forcing. The 
poloidally elongated structures can experience local forcing in the radial direction. Diagonal 
connections in the {vy){x,t) spatiotemporal GAM pattern appear through local radial zonal 
jet migration if the life time of the jets is comparable to the GAM period. For lower cb 
(e.g. 0.1) the vortex turnover period is much smaller than the GAM modulation time of the 
zonal flows, and no diagonal motion is observed in spatiotemporal diagrams. 

Up-down asymmetry will not change the average effect of geodesic forcing on [ky = 0) 
zonal flows and GAMs. It may however have an indirect effect on poloidally localized 
[ky <^ 1) zonal jets. Asymmetric shaping also influences local magnetic shear, which is, 
for example, strongly enhanced near a (lower) single-null X-point [36]. Zonal jets emerging 
through straining of drift wave vortices can be weaker at the (downward) poloidal location 
of strong magnetic shear, and may experience less geodesic forcing than on the other (up) 
side. Under such circumstances an effective geodesic forcing of zonal jets might actually 
appear. In order to consistently account for local magnetic shear effects in interplay with 
zonal flow shearing a 3-D (flux-tube) model is indispensible. This case is therefore not 
further considered within the simplified 2-D model of the present discussion. 
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It has been demonstrated in the preceeding sections that unphysical radial propagation 
of zonal flows and GAMs by geodesic forcing may apparently occur for certain computa- 
tional set-ups. Radial propagation of zonal modes by other possible causes (for example by 
inhomogeneous FLR effects) is not disputed here. 

IX. "SECOND SOUND" IN DRIFT WAVE TURBULENCE 

Radially propagating zonal modes, however they may originate, bear a striking resem- 
blence to the phenomenon of "second sound" in quantum codensates. Second sound in 
two-fluid systems, for example in superfluid He, describes the transport of temperature and 
entropy in the form of waves rather than by diffusion [44, 45]. Phonons ("first sound") are 
modulated by temperature waves ("second sound"), while the normal component and the 
condensate oscillate out of phase. 

Drift wave turbulence can be considered within the framework of wave kinetics [46] . Wave 
density and temperature in this context do not refer to the fluid moments of the plasma 
particles, but rather to the population density and wave energy of the drift mode excitations. 

The spectral drift wave mode density of the turbulence at a specific radial location is 
modulated in time by the passing wave of the radially propagating zonal flow condensate. 
Although the condensate phase docs not by itself participate in radial plasma transport this 
leads to a radial wave like transfer of spectral mode energy, when the turbulent phase is 
co-existing with the propagating condensate. Phenomenologically this propagation shows 
characteristics of a classical analogon to second sound in quantum condensates: in a wave 
kinetic picture a mode temperature of the turbulent phase may be defined by the mean 
spectral energy, which is modulated by zonal condensate mode propagation in the form of 
a radial wave rather than by diffusive spreading. 

The source of drift wave turbulence is, with respect to its two-dimensional character 
and dual cascade property, located within the spectral distribution around k±ps ~ 1 and 
is for many cases well separated from both dissipation and modulation (condensate) scales. 
Weakly damped second sound in drift wave turbulence of this form appears possible, in 
contrast to general 3-D turbulent Kolmogorov type systems [47]. A future more rigorous 
theoretical consideration of this phenomenon surely is eligible. 
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X. CONCLUSIONS 



Basic mechanisms of local and global radial geodesic forcing on zonal modes have been 
discussed with the 2-D HW standard model of drift wave turbulence. The main results are 
summarized as follows: 

• Radial propagation oi ky <^ 1 zonal modes can occur by local geodesic forcing, which 
has been demonstrated by poloidally local simulations. 

• A net radial geodesic forcing oi ky — zonal modes (zonal flows and GAMs) can not 
be expected in any closed toroidal magnetic confinement configuration, as the flux 
surface average of geodesic curvature always identically vanishes. 

• Unphysical radial propagation of zonal flows and GAMs may appear in simulations 
that use inappropriately simple geometry models violating this zero average property. 

• Apparent radial motion of spatio-temporal zonal mode patterns may be visually diag- 
nosed in the presence of geodesic acoustic oscillations. 

• Zonal jets with ky <^ 1 occur as vortices strained out by ky — zonal flows or GAMs 
and may experience local geodesic forcing. 

• Radially migrating poloidally localized zonal jets lead to a diagonal connection in 
spatio-temporal GAM patterns, and in conclusion to an illusion of radial GAM prop- 
agation. 

• Up-down asymmetry does not alter the average geodesic curvature from zero, and 
therefore does not affect radial geodesic forcing of zonal flows and GAMs. 

• Up-down asymmetry may influence the local amplitude of zonal jets through local 
magnetic shearing. 

• Phenomenologically a radial propagation of zonal modes shows some characteristics 
of a classical analogon to second sound in quantum condensates. 

While the applicability of the 2-D HW model on realistic fusion plasmas is definitely 
limited, it has served to elucidate basic physical effects of the magnetic field geometry on 
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drift wave turbulence and zonal flows. Any future computational studies with more complete 
models on radial propagation of zonal modes should not rely on inappropriately simplified 
geometry models but have to make use of actual MHD equihbria. 
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